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I.  INTRODUCTION 


A.  WHAT  IS  SPUTTERING? 

The  ejection  of  surface  atoms  from  solid  surfaces  under 
bombardment  by  energetic  particles  is  known  as  sputtering. 
The  particles  may  be  ions,  neutral  atoms,  neutrons, 
electrons  or  photons.  Ejection  is  caused  by  collisions 
between  the  incoming  particles  and  atoms  in  the  selvage, 
i.e.,  the  surface  layers  of  the  solid.  The  processes 
involved  in  sputtering  are,  in  principle,  similar  to  those 
causing  radiation  damage  in  the  bulk  of  a  solid.  These 
processes  occur  far  from  thermal  equilibrium,  which  means 
that  sputtering  is  different  from  evaporation,  which  does 
occur  at  thermal  equilibrium. 

Generally,  an  incoming  particle  will  collide  with  the 
atoms  of  the  solid,  thereby  transferring  energy  to  the 
atomic  nuclei.  If  the  energy  transferred  to  an  atomic 
nucleus  is  greater  than  the  binding  energy  at  the  lattice 
site,  a  primary  knock-on  atom  (PKA)  is  created.  The  PKA 
will  collide  with  other  target  atoms  distributing  the 
energy  by  way  of  a  collision  cascade.  An  atom  is  sputtered 
if  the  energy  transferred  to  it  has  a  component  normal  to 
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the  surface  which  is 


larger  than  the  surface  potential 


energy  barrier1 . 

One  measure  of  sputtering  is  the  sputtering  yield,  Y, 
defined  as  the  mean  number  of  atoms  removed  per  incident 
particle.  In  counting  the  removed  atoms,  only  those  of  the 
solid  are  included,  while  reflected  ions  or  incident  atoms 
that  have  been  reemitted  are  not  taken  into  account.  The 
particles  removed  by  sputtering  are  emitted  with  a  broad 
distribution  of  energy,  in  different  excitation  and  charge 
states,  and  at  all  exit  angles.  Differential  yields  are 
defined  to  describe  the  variation  of  the  yield  with  these 
parameters. 

B.  APPLICATIONS  OF  SPUTTERING 

Sputtering  was  long  regarded  as  an  undesirable  effect 
which:  i)  destroys  cathodes  and  grids  in  gas  discharge 
tubes  or  ion  sources,  ii)  contaminates  a  plasma  and 
surrounding  walls,  making  it  a  major  research  area  in 
fusion  technology,  and  iii)  causes  the  destruction  of 
diaphragms  and  targets  in  accelerators  and  in  high-voltage 
electron  microscopes.  However,  sputtering  is  used  today 
for  many  applications  and  has  become  an  indispensible 
process  in  modern  technology. 

1  The  potential  energy  of  a  surface  layer  atom  is 
normally  referred  to  as  the  surface  binding  energy  (SBE). 
The  SBE,  however,  only  has  meaning  in  an  undamaged  lattice 
and  does  not  define  the  energy  required  to  remove  a  surface 
atcm  once  the  collision  process  has  begun. 
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Sputtering  allows  a  controlled  removal  of  tightly  bound 
surface  layers  on  a  nearly  atomic  scale  and  possible 
submicron  spacial  resolution  if  a  well  focused  or  rastered 
beam  is  used-  The  atoms  removed  can  be  analyzed  in  a  mass 
spectrometer,  which  gives  information  about  surface 
concentration,  or  can  provide  a  depth  profile  if  the 
surface  is  continuously  bombarded. 

Both  the  removal  of  atoms  from  a  surface  and  the  flux 
of  atoms  leaving  the  surface  are  successfully  applied,  for 
example:  in  sputter  ion  sources,  to  obtain  atomically  clean 
surfaces;  in  micromachining;  and  for  depth  profiling  of 
thin  films  and  surfaces-  One  of  the  largest  applications 
of  sputtering  is  the  deposition  of  thin  films  on  a  large 
variety  of  substrates.  These  films  may  be  several  square 
meters  in  area,  or  may  be  extremely  small  as  in 
microelectronics. 

C-  HISTORICAL  OVERVIEW 

Sputtering  has  been  investigated  in  most  detail  for 
bombardment  of  monatomic  solids  with  mercury  ions,  with 
noble  gas  ions  and  with  hydrogen  ions.  Sputtering  yields 
may  lie  between  0  and  104  atoms  per  incident  ion,  but  are 
typically  one  to  five.  The  yields  depend  on:  i)  the 
incident  particle's  energy,  mass  and  angle  of  incidence; 
ii )  the  mass  of  the  target  atoms;  iii)  the  range  of  crystal 
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order6  and  the  crystal  orientation  of  the  solid  and  iv)  the 
surface  binding  energies  of  the  target.  The  yield  is, 
however,  nearly  independent  of  the  temperature  [Refs.  1,  2, 
3].  Below  a  threshold  energy,  which  is  about  20-40  eV  for 
normally  incident  ions,  no  sputtering  takes  place.  Above 
this  threshold  the  yield  increases  with  incident  ion  energy 
and  reachs  a  broad  maximum  in  the  energy  region  of  5-50 
keV.  At  higher  energies  the  sputtering  yield  decreases. 
This  decrease  is  related  to  the  larger  penetration  of  the 
ions  into  the  solid  and  the  lower  energy  deposition  into 
the  surface  layers. 

Several  versions  of  the  theory  of  sputtering  currently 
exist  in  the  literature,  each  of  which  exhibits  the 
propensities  of  their  authors.  Each  version  has  some 
experimental  support,  but  certain  approaches  are  in  better 
repute  than  others.  This  review  makes  no  attempt  to 
rationalize  the  various  approaches,  but  rather  it  is 
concerned  with  providing  a  historical  overview  of  the 
development  of  the  theories. 

Sputtering  was  first  recorded  in  1852  by  Grove  [Refs. 
4,  5]  who  noticed  the  disintegration  of  cathodes  in  glow- 
discharge  tubes  with  subsequent  deposition  of  the  material 
on  the  surrounding  glass  walls.  This  became  known  as 


6  The  range  of  crystal  order  refers  to  the  distance 
over  which  the  crystal  maintains  the  regular,  repetetive 
gridlike  arrangement  of  atoms,  known  as  a  lattice. 
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"cathode  sputtering"  in  the  literature.  It  took  almost 
fifty  years  before  Goldstein  [Ref. 63  presented  compelling 
evidence  that  the  sputtering  effect  was  caused  by  positive 
ions  of  the  discharge  striking  the  cathode. 

Unfortunately,  early  investigations  of  sputtering  were 
wasted  before  it  was  realized  that  the  pressure  surrounding 
the  target  was  a  critical  parameter  and  had  to  be  low 
enough  to  allow  the  sputtered  atoms  to  escape  the  surface. 
Many  anomalies  in  early  experimental  results  can  be 
accounted  for  by  the  simple  fact  that  the  mean  free  path  of 
the  sputtered  particles  was  less  than  the  distance  to  some 
"collector  surface",  which  allowed  the  liberated  particles 
to  diffuse  back  to  the  surface  after  collision  with  a  gas 
atom.  Penning  and  Moubis  [Ref.  73  were  the  first  to 
conclusively  demonstrate  the  effect  of  ambient  pressure, 
even  though  the  view  that  pressure  was  an  important  factor 
had  been  held  for  a  few  years  prior  to  their  publication. 

The  work  of  Gunther schultze  and  Meyer  [Refs.  8,  93  is 
one  of  the  few  early  experiments  preceding  the  findings  of 
Penning  and  Moubis  which  satisfied  the  conditions  for  a 
reproducible  sputtering  determination. -  Gunther schul tze  and 
Meyer  operated  their  discharge  tube  at  relatively  high 
vacuum  (for  that  era)  with  a  copper  target  and  took  the 
precautions  of  removing  initial  layers  of  the  target 
material. 
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The  elimination  of  adsorbed  gases  and  surface  oxide 
layers  as  a  critical  factor  in  determining  the  true  yield 
was  further  described  by  Arifov,  et  al.  ,  in  1963  [Ref.  10] • 
Yonts  and  Harrison  also  presented  evidence  [Ref.  11]  that 
surface  recont amination  from  background  gases  was  a 
significant  factor  in  quantitative  sputtering  yield 
measurements.  Criteria  were  therefore  developed  for 
conducting  "clean”  experiments.  For  a  single  crystal  metal 
this  means  the  removal  of  surface  and  near-surface  damage 
and  impurities  by  mechanical  abrasion  and  polishing,  with 
thermal  annealing  and  "sputter  cleaning"  used  to  remove  any 
residual  damage.  It  is  possible  to  go  to  greater  lengths 
to  ensure  damage-  and  cont aminan t - f r ee  surfaces,  for 
example,  using  low  energy  heavy  ion  bombardment  to  remove 
surface  oxide  (conducted  in  ultrahigh  vacuum)  followed  by 
thermal  annealing  to  remove  surface  damage  caused  by  the 
ion  bombardment.  Information  about  the  atomic  structure  of 
the  surface  can  then  be  obtained  by  low  energy  electron 
diffraction  ( LEED ) ,  field  ion  microscopy  of  the  new 
tunneling  electron  microscopy.  Subsequent  sputtering 
experiments  could  then  be  conducted  in  an  ultrahigh  vacuum 
environment  (less  than  10-t0  Torr)  to  preclude  adsorption 
of  gas  molecules  on  the  surface. 

The  concept  of  an  individual  sputtering  event  on  an 
atomic  scale  was  first  recognized  and  analyzed  by  Stark. 
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He  originated  the  so-called  hot-spot  model3  of  sputtering 
[Ref.  123,  and  subsequently,  a  collision  theory  viewing 
sputtering  as  a  sequence  of  binary  collision  events 
initiated  by  one  bombarding  ion  at  a  time  [Ref.  131. 
Utilising  the  conservation  laws  of  the  theory  of  elastic 
collisions  and  the  concept  of  collision  cross  sections  in 
his  analysis,  Stark  correctly  interpreted  the  observed 
energy  dependence  of  the  sputtering  yield  of  hydrogen  ions 
bombarding  a  metallic  target. 

Subsequent  investigators  took  a  different  stand  on  the 
issue.  After  Kingdon  and  Langmuir  successfully  applied 
Stark's  collision  theory  [Ref.  143,  the  erroneous 
impression  arose  that  the  collision  theory  would  imply  that 
sputtering  was  a  single-collision  process  resulting  in  a 
strongly  peaked  angular  distribution  of  the  flux  of 
sputtered  particles.  Only  after  the  demonstration  by 
Wehner  [Refs.  15,  lb]  of  crystal  structure  effects  in  the 
flux  of  sputtered  particles  did  it  become  evident  that 
local  evaporation  alone  could  hardly  explain  the  sputtering 
phenomena.  Wehner  discovered  that  the  angular 
distributions  of  the  atoms  ejected  from  single  crystal 
targets  showed  maximum  intensities  in  directions 


3  The  hot-spot  model  was  considered  to  be  the 
evaporation  of  target  material  from  a  microscopically  small 
region  with  a  high  local  temperature  due  to  individual  ion 
impact. 
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corresponding  to  the  more  closely  spaced  atomic  rows 
passing  through  the  target  surfaces. 

The  discovery  of  such  ejection  patterns  or  "Wehner 
spots0  revived  the  interest  in  the  collision  theory  of 
sputtering.  Independent ly ,  Keywell  [Refs.  17,  183  made  a 
first  attempt  to  formulate  Stark's  multiple-collision  model 
in  terms  familiar  from  neutron  transport  theory.  Keywell's 
work,  as  well  as  subsequent  calculations  by  Harrison  [Ref. 
193  were  important  steps  in  the  sense  that  probability 
concepts  as  expressed  by  collision  cross  sections  finally 
made  their  entrance  into  sputtering  theory. 

The  theory  of  collision  cascades  has  been  further 
advanced  by  the  work  of  Leibfried  [Ref.  20],  Lindhard,  et 
al.  CRef.  21],  Dederichs  [Ref.  22],  Robinson  [Ref.  23], 
Thompson  [Ref.  24],  and  Sigmund  [Refs.  25,  2b],  who  has 
also  collected  the  available  knowledge  in  a  transport 
theory  of  sputtering  [Ref.  27].  In  an  amorphous  solid  the 
collision  cascades  have  been  described  by  a  Boltzmann 
transport  equation.  Sigmund  was  able  to  obtain  first  order 
asymptotic  solutions  for  the  linear  cascade  regime4 . 
However,  only  polycrystalline  solids  with  randomly  oriented 
crystallites  can  be  approximated  by  an  amorphous  solid. 


4  The  linear  cascade  regime  applies  to  bombardment 
with  ions  of  medium  to  large  atomic  number  in  the  keV 
energy  region.  Here  large  collision  cascades  can  develop, 
however,  moving  target  atoms  are  considered  to  collide  only 
with  target  atoms  at  rest. 
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Despite  its  shortcomings,  this  theory  has  become  successful 
as  a  reference  standard  for  sputtering  yield  measurements. 
It  has  also  been  a  reasonable  starting  point  for  further 
refinements  of  the  theory.  During  the  rapid  development  of 
the  theory  of  collision  cascades,  the  hot-spot  (or  "spike”) 
theory  has  received  little  attention  despite  repeated 
claims  of  a  thermal  component  observed  in  energy  spectra 
of  sputtered  atoms  [Ref.  28]. 

In  the  period  from  1955  to  1965,  the  dominant  effort  in 
the  study  of  sputtering  on  metals  was  actually  spent  in  the 
investigation  of  crystal  lattice  effects.  These  effects 
fell  into  two  distinct  groups,  the  channeling  theory  of  the 
sputtering  yield  and  the  focusing  collision  theory  of 
ejection  patterns.  The  channeling  theory  accounts  for  the 
angular  variations  of  the  yield  cloee  to  the  principle 
crystal  axes  and  planes,  but  is  not  very  successful  for 
more  general  directions.  The  focusing  collision  theory 
gives  a  qualitative  description  of  the  ejection  pattern 
spots  from  cubic  metals  at  high  bombarding  energies,  but 
cannot  successfully  interpret  low  energy  irradiation 
results  or>  observations  on  hexagonal  or  other  low  symmetry 
materials.  The  discovery  of  spot  patterns  by  Wehner  [Refs. 
15,  16]  inspired  Silsbee's  iden t if icat ion  [Ref.  29]  of  the 
focusing  collision  sequence  as  a  means  for  long-range 
transport  of  momentum  in  crystals  at  low  energy.  An 
alternative  to  the  focusing  theory  of  the  ejection  pattern 
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spate  was  proposed  by  Lehmann  and  Sigmund  [Ref.  301.  Their 
model  stipulates  that  the  target  surface  have  an  ordered 
structure,  but  not  necessarily  long  straight  rows  of  atoms 
intersecting  the  surface.  This  model  emphasizes  the  role 
of  surface  structure  and  of  the  surface  binding  energy  in 
monocrystal  sputtering. 

D.  COMPUTER  MODELING 

In  the  past,  sputtering  theory  needed  only  to  predict 
the  sputtering  yield  as  a  function  of  energy  for  a  given 
incident  particle  and  polycrystalline  target  of  a  single 
metal.  The  experimental  environment  led  the  theorists  to 
statistical  theories  involving  ensemble  averages  over 
distribution  functions,  ultimately  culminating  in  the  work 
of  Thompson  [Ref.  24]  and  Sigmund  [Ref.  25]  which  predict 
the  yield  reasonably  well.  However,  the  chief  limitation 
of  these  theories  is  their  inability  to  deal  with 
experimental  results  attributable  to  the  regular  structure 
in  a  single  crystal  target.  As  is  the  case  with  most 
theory,  the  statistical  sputtering  theories  are  also 
constrained  by  the  lack  of  analytic  tract ability . 

With  the  advent  of  high  speed  computing  capabilities  it 
became  practical  to  include  the  full  three-dimensional 
structures  of  crystals  in  the  theoretical  models  of 


sputtering. 

The 

investigations 

of 

Vineyard 

and  his 

colleagues 

[Ref.  31] 

demonstrated 

the 

power  of 

computer 

simulations  in  isolating  and  under st anding  the  elementary 
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processes  of  defect  production  and  migration  in  radiation 
damage.  Such  techniques  were  first  applied  to  monocr y st al 
sputtering  problems  by  Harrison  and  his  co-workers  [Ref. 
32].  They  proposed  ejection  mechanisms  to  account  for  the 
Wehner  spots  which  differ  from  both  the  focusing  collision 
theory  and  the  model  of  Lehmann  and  Sigmund. 

These  computer  simulations  can  be  divided  into  two 
categories*  Historically,  the  binary  collision  ( BC ) 
simulation  models  were  developed  first*  The  binary 
collision  simulation  is  known  as  an  "event  store”  model. 
The  program  follows  a  cascade  of  moving  particles 
sequentially,  moving  only  one  particle  at  a  time,  but 
remembering  events,  i. e. ,  collisions,  as  it  progresses. 
The  models  usually  invoke  the  binary  collision 
approximation  [Ref.  35J.  This  assumes  that  each  particle 
interacts  with  only  one  other  particle  at  any  ti me,  with 
the  other  particle  usually  considered  to  be  stationary. 
Thus  these  models  are  inherently  linear  calculations  and 
therefore  they  have  many  of  the  same  restrictions  as  the 
linear  transport  calculations. 

Initially,  these  models  simply  mimicked  the  statistical 
theories,  then  Robinson  and  Oen  included  the  target  crystal 
structure  [Ref.  331  and  effectively  discovered  channeling. 
The  channeling  concept  fed  back  into  sputtering  theory  and 
Onterder linden  [Ref.  34]  introduced  the  "channeling" 
sputtering  model,  which  satisfactorily  explained  the 
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relative  magnitude  of  the  sputtering  yield  from  single 
crystal  target  surfaces  at  low  keV  ion  energies.  Since 
their  inception,  the  BC  programs,  MARLOWE  [Ref.  35],  TRIM 
[Ref.  36]  and  MORLAY  [Ref.  37],  have  assumed  that  the 
physical  model  was  well  understood  and  that  the  simulation 
should  be  designed  to  compute  results  from  it  efficiently. 
Each  program  has  its  particular  strengths  and  efficiencies, 
because  each  was  designed  to  address  specific  criteria. 

The  second  method  of  simulation  is  the  multiple 
interaction  (Ml)  simulation  and  is  based  on  molecular 
dynamics.  MI  simulation  was  motivated  by  the  work  uf 
Gibson,  Goland,  Milgram  and  Vineyard  [Ref.  31],  with  the 
proponents  of  this  model  having  concentrated  on 
understanding  the  fundamentals  of  ejection  processes. 

This  type  of  simulation  is  known  as  an  "timestep"  model 
and  involves  solving  Newton's  laws  of  motion,  often 
expressed  in  Hamiltonian  form,  for  many  particles.  These 
calculations  proceed  sequentially  in  time,  following  the 
collisions  of  many  particles  simultaneously.  An  advantage 
of  the  model  is  that  there  are  no  built  in  constraints  as 
to  the  behavior  of  a  collision  cascade  and,  given  enough 
cpu  time  and  memory,  will  accurately  follow  an  ion  induced 
cascade  to  its  ultimate  ther malizat ion. 

A  new  development  is  that  of  a  hybrid  between  the  two 
types  of  codes  by  Harrison,  et.  al.  [Ref.  38].  The  goal  of 
which  is  to  have  a  simulation  which  retains  the  "exact” 
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treatment  of  the  time  step  model  in  the  surface  region  and 
minimizes  computation  time  for  the  remainder  of  the 
cascade . 


As  the  simulations  have  become  increasi 
sophisticated;  attention  has  turned  to 


semiconductor  materials 
distinct  problem  due 
potentials.  Computer 
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3  Private  communica t ions  between  Professor  B.  J. 
Garrison  of  Pennsylvania  State  University  and  Professor 
D.  E.  Harrison,  Jr.  of  the  Naval  Postgraduate  School. 
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II 


OBJECTIVES 


Numerous  experimental  and  theoretical  investigations  of 
the  energy  and  angular  distributions  of  atoms  and  molecules 
ejected  from  single  crystal  surfaces  due  to  ion  bombardment 
have  been  made  over  the  years.  Wehner  was  the  first  to 
observe  spots,  or  anisotropic  angular  distributions  [Ref. 
15,  163.  His  method  of  detection,  however,  did  not  allow 
him  to  distinguish  the  mass  or  energy  of  the  ejected 
species.  In  more  recent  times  energy  and  angular 
distributions  of  ions  have  been  measured  by  angle-resolved 
Secondary  Ion  Mass  Spectroscopy  (SIMS)  [Refs.  39,  40,  41, 
42].  The  theoretical  investigations,  mainly  computer 
simulations,  have  been  used  to  predict  the  distribution  of 
neutral  species  [Refs.  43,  44,  45,  46].  These  molecular 
dynamics  calculations  of  the  ion  impact  event  have  helped 
to  i)  formulate  the  mechanisms  of  particle  emission  from 
single  crystal  surfaces,  ii )  provide  an  atomistic  picture 
of  the  damage  created  within  the  material  itself,  and  iii) 
unravel  the  origin  of  angular  anisotropies. 

Experimental  measurements  of  particle  trajectories  have 
not  been  possible  at  the  level  necessary  to  make  direct 
comparisons  to  theoretical  predictions.  However, 
qualitatively  and  semi -quant it at ively  favorable  comparisons 
have  been  made  between  the  predicted  neutral  distributions 
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and  the  measured  ion 


42] 


distributions  [Refs.  40,  41, 

Rigorous  quantitative  comparisons  are  difficult  since  the 
ionization  probability  is  not  known.  The  ionization 

probability  is  required  since  the  ions  experience  a 
different  interaction  potential  than  the  neutrals  as  they 
leave  the  surface.  In  particular,  it  would  be  advantageous 
to  be  able  to  directly  compare  experimental  and  theoretical 
distributions  in  order  to  use  the  measured  distributions  to 
gain  structural  information  about  the  surface  and  to 
improve  the  theoretical  model  so  as  to  gain  a  better 
under standing  of  the  scattering  process. 

With  the  recent  advent  of  the  energy-  and  angle- 
resolved  neutral-particle  (EARN)  detector  by  N.  Winograd, 
et.  al.  [Ref.  47],  distributions  for  rhodium  atoms  ejected 
from  clean  and  oxygen  covered  Rh(lll}  surfaces  due  to  ion 
bombardment  have  been  obtained  [Refs.  43,  49,  50].  Energy 

and  angular  distributions  reported  by  this  method 
constitute  the  first  reported  neutral -particle  ejection 
distributions  for  material  desorbing  from  genuinely  clean 
well-defined  surfaces.  Surface  cleanliness  and  structure 
were  monitored  by  low  energy  electron  diffraction  ( LEED ) 
and  by  Auger  electron  spectroscopy  (AES)  as  well  as  by  the 
reproducibility  of  the  EARN  measurements. 

The  results  of  computer  simulations  of  argon  ion 
bombardment  of  RhClll)  still  do  not  correlate 

quantitatively  with  the  reported  EARN  results.  This  lack 
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experimental  results,  this  being  left  for  future 
when  more  experimental  data  becomes  available. 


studies 
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III. 


THE  MODEL 


A.  MULTIPLE  INTERACTION  (Ml)  SIMULATION 

The  MI  simulation  used  to  model  the  ion  bombardment 
events  is  QDYN,  developed  by  Professor  Don  E.  Harrison,  Jr. 
at  the  Naval  Postgraduate  School.  A  detailed  description 
of  the  computer  simulation  method,  the  mathematical  model 
it  is  based  on,  and  the  integration  scheme  have  been 
published  previously  by  Harrison,  et.  a 1 .  [Refs.  32,  53, 
54  3  . 

Briefly,  QDYN  is  a  sputtering  simulation  which  proceeds 
through  a  series  of  timesteps,  each  of  which  consists  of 
the  following:  i)  summation  of  the  pairwise  forces  for  each 
atom;  ii)  calculation  of  new  velocities  and  positions  at 
the  end  of  the  timestep;  iii)  movement  of  atoms  to  their 
new  positions;  and  iv)  test  of  energy  conservation. 

QDYN  computes  the  interaction  forces  by  performing  a 
numerical  solution  of  Hamilton's  classical  equations  of 
motion,  using  a  predictor -cor rector  integration  scheme 
which  bases  the  time  increment  of  each  step  on  the  fastest 
moving  particle.  This  differs  from  the  form  of  most  other 
numerical  integrators  because  position  and  velocity  are 
calculated  for  the  same  instant  of  time.  More 
sophisticated  integration  schemes  do  not  improve  the 
calculations  and  only  increase  the  computation  costs . 
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The  system 

consists 

of  ah 

incoming  argon 

ion  and  a 

single-crystal 

rhodium 

target 

containing  up  to 

1854  atoms. 

Each  particle  is  char act er ized  by  i)  its  mass,  ii)  the 


force  laws  by 

which  it 

interacts  w 

ith  other  particles, 

iii  ) 

its  position, 

and  iv) 

its 

velocity 

Once  the  incoming 

ion 

strikes  the 

target 

at 

a  spec 

ified  impact  point 

the 

positions  and 

momenta 

of 

the  ion  a 

nd  all  the  atoms  in 

the 

single-crystal 

target 

are 

allowed 

to  develop  in  time. 

This 

calculation,  called  a 

trajectory, 

is  a  simulation 

of  a 

single  collision  cascade.  The  collision  cascade  is 
terminated  when  the  momentum  has  dissipated  through  the 
lattice  and  no  more  atoms  can  be  ejected.  In  this  metallic 
rhodium  system,  this  corresponds  to  a  maximum  particle 
energy,  in  the  solid,  of  2. 0  eV.  Finite  lattice 
temperature  effects  and  energy  losses  due  to  electronic 
excitations  have  not  been  included  in  this  calculation. 

Two  ancillary  programs  were  used  to  investigate  each 
simulation.  ANMQL  analyzes  the  final  positions  and  momenta 
to  determine  yields,  the  energy  distribution,  possible 
multimer  formation,  and  the  angular  distribution  of  the 
ejected  atoms,  all  of  which  are  experimentally  measurable 
quantities.  In  addition,  information  about  atoms  ejected 
per  single  ion  (ASI)  and  ejection  times  is  provided. 
Graphical  representations  of  the  distributions  and 
quantities  mentioned  above  were  provided  by  ANPLOT.  An 
alternate  output  available  from  GDYN  allows  sufficient 
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Information  to  be  printed  out  so  that  the  collision  cascade 


can  be  traced. 

B.  PHYSICAL  CONSIDERATIONS 

1 .  Lattice  Orientation  and  Constants 

A  Rh {111)  orientation  was  chosen  to  match  the 
conditions  of  the  energy-  and  angle-resolved  neutral- 
particle  (EARN)  experiments  [Refs.  48,  49,  50].  Rhodium  is 

a  face  centered  cubic  (fee)  metal  with  a  lattice  constant 
(a0)  of  3. 804A.  Other  physical  constants  of  rhodium  used  in 
the  development  of  a  Rh-Rh  potential  function,  include:  the 
bulk  heat  of  atomization  (5.76  eV )  and  the  bulk 

compressibility  (3.  626x10**  7  cm2  /  kg). 

The  basic  distance  unit  defined  for  this  simulation 
is  the  lattice  unit  ( LU ) ,  which  is  one  half  the  lattice 
constant  for  cubic  materials. 

1  LU  =  0. 5  a*  =  1. 902A 

2.  Lattice  Size  and  Containment  Criteria 

To  model  the  5  keV  bombardment  of  the  EARN 
experiments  and  describe  most  of  the  ejection  events  would 
require  a  target  size  that  is  computationally  not  feasible. 
In  order  to  gain  the  most  complete  information  about  the* 
Rh ( 1 1 1 } / Ar +  system  this  study  conducted  its  computations  at 
ion  bombardment  energies  of  500  eV  and  2  keV.  These 
energies  and  their  corresponding  lattice  size  were  chosen 
in  order  to  i)  optimize  containment  of  the  collision 
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cascade  within  the  lattice,  ii )  maximize  the  number  of 
atoms  sputtered  to  ensure  statistical  confidence  in  the 
results,  and  iii )  minimize  the  computational  time  necessary 
to  complete  a  simulation. 

The  inability  of  the  lattice  of  target  atoms  to 
contain  the  entire  cascade  created  by  the  incident  ion  is 
called  a  failure  of  "containment" .  In  principle,  there  is 
no  problem  for  a  moderate-size  computer  to  store  the 
information  for  even  tens  of  thousands  of  atoms,  but,  in 
practice,  running  time  considerations  (i.e.,  money)  limit 
the  target  to  roughly  2000  atoms.  Cascades  created  by  keV 
ions  are  known  to  involve  many  more  than  2000  atoms,  so  the 
problem  seems  insoluble.  After  long  experience  with  these 
computat ions,  Harrison  has  defined  a  more  practical 
definition  of  containment  [Ref.  55]:  "If  an  increase  in 
target  size  does  not  change  the  results  of  the  computation, 
the  trajectory  is  effectively  contained  for  the  purposes  of 
that  computation". 

In  practice,  for  any  observable,  one  finds  that 
containment  is  approached  asymptotically  as  the  target  size 
increases.  Harrison  also  states  that: 


As  one  might  anticipate  from  experimental  experience 
with  absolute  values,  yield  provides  the  most  severe 
test  of  containment;  so  absolute  yield  containment 
guarantees  everything  else.  As  the  target  size 
increases,  all  other  global  results  reach  constant 
values  for  much  smaller  targets  than  those  required  to 
produce  stable  yield  values. 
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This  fact  is  born  out  by 


Figure  1,  which  is  a 


graphical  summary  of  numerous  calculations  done  on  clean 
Rh(lll>  at  various  ion  energies  and  lattice  sizes.  These 
calculations  provide  the  baseline  information  on  clean 
Rh(lll}  as  well  as  acting  as  verification  of  containment. 

The  following  lattice  sizes  (see  Table  1) 
correspond  to  the  information  provided  in  Figure  1.  The 
lattice  dimensions  are  given  in  terms  of  the  number  of 
lattice  planes  in  the  x,  y,  and  z  directions.  Figure  2 
provides  a  graphical  description  of  the  lattice  plane 
numbering  system. 


TABLE  1 

LATTICE  DIMENSIONS 


Lat  t  i 


:e 

Di 

mension 

#  Rh  Atoms  in  the  Lattice 

18 

X 

5 

X 

12 

541 

ia 

X 

6 

X 

12 

848 

21 

X 

5 

X 

13 

883 

21 

X 

6 

X 

13 

820 

23 

X 

a 

X 

15 

1038 

25 

X 

6 

X 

17 

1278 

27 

X 

a 

X 

19 

1540 

25 

X 

a 

X 

17 

1701 

The  lattice  sizes  chosen  for  the  computations  were: 

Ion  Energy  Lattice  Size 

500  eV  18  x  5  x  12 

2000  eV  25  x  6  x  17 

Even  though  appreciable  movement  of  the  edge  atoms 

was  observed  prior  to  the  last  atom  leaving  the  surface, 

the  significant  features  of  the  dynamics  were  not  altered 
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when  a  larger  lattice  size  was  used.  All  simulations  were 
run  on  an  IBM  3033  at  the  Naval  Postgraduate  School  with 
computation  times  on  the  order  of  1  CPU  hour  for  the  500  eV 
simulations,  and  3  CPU  hours  for  the  2  keV  simulations. 

For  the  13x5x12  lattice  there  are  103  atoms  in  the 


first 

layer  and 

for 

the  25x6x17 

lattice  there  are 

213 

atoms . 

To  aid  in 

de 

scribing  the 

simulation  events 

all 

atoms  < 

and  ions) 

in 

the  simulation 

are  numbered,  with 

the 

ion  bei 

ng  number  1 

and 

the  lattice  atoms  numbered  beginn 

i  ng 

with  number  2  at  the  origin  and  progressing  upward  across 
the  rows,  one  layer  at  a  time.  Figures  3  and  4  delineate 
the  numbering  system  for  the  first  layer  of  each  lattice 
system. 

Another  containment  considers tion  deals  with  those 
atoms  that  leave  the  micr ocryst al 1 i te  through  the  bottom 
and  sides.  The  simulation  program  considers  these  atoms  to 
be  lost  to  the  collision  process,  but  continues  to  track 
them  for  energy  conservation  purposes  until  the  model  is 
terminated.  This  is  justified  for  those  atoms  that  depart 
the  lattice  through  the  bottom,  since  almost  all  of  the 
collision  events  that  contribute  to  sputtering  occur  in  the 
first  four  layers  of  the  lattice.  The  atoms  that  are 
ejected  from  the  sides  of  the  crystal  are  ignored  in  order 
to  improve  computational  efficiency.  Harrison,  through 
long  experience  with  this  model  has  determined  that  the 
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deletion  of  these  atoms  has  little  impact  on  the  number  of 


atoms  sputtered,  altering  the  results  by 
Alternate  approaches  to  those  a 
the  sides  of  the  microcrystal  would  be: 
at  the  side  boundaries,  or  ii)  to  reintr 
opposite  side  of  the  lattice  from  which 
the  same  momenta. 

3.  Statistical  Ensemble 

An  ensemble  average  is  an  averag 
over  a  group  of  similar  systems.  S 
ensemble  averages,  where  each  trajectory 
an  ensemble  of  all  possible  trajectories 
particular  atom  ejection.  For  a  monatomi 
point  on  the  target  surface  leads 

trajectory.  Using  point  and  line  symmetr 
one  establishes  a  n? pr e s e n t a  t  i  v-t?  a r tea  i n 
is  an  elementary  symmetry  zone  of  the 
lattice  orientation.  An  ensemble  of 

developed  by  uniform  sampling  of  the  repr 
For  our  purposes  a  rectangular  re 
of  dimensions  0.707  LU  by  1.225  LU,  wa 
origin  placed  at  the  cegter  of  a  target 
4).  Uniform  sampling  of  the  re pres 

accomplished  by  designating  300  impact 
rectangle  (Fig.  5). 


about  0. IX. 
toms  lost  through 
i )  to  reflect  them 
oduce  them  on  the 
they  exited,  with 

e  of  any  property 
irnulations  produce 
is  a  sample  from 
realizable  from  a 
c  ion,  each  impact 
to  a  different 
y  transformations, 
the  target,  which 
surface  of  that 
trajectories  is 
esentative  area, 
present ative  area, 
s  chosen,  with  its 
atom  (Figs.  3  and 
entative  area  was 
points  within  the 
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In  order  to  reduce  computation  time,  the  necessity 
of  using  a  uniform  distribution  of  300  trajectories  was 
investigated.  The  impact  points  were  divided  into  two 

groups  of  150  points  each  (indicated  by  x's  and  +■  '  s  in 
Fig.  5)  and  the  results  of  separate  simulations  of  500  eV 
Ar+  on  clean  Rh(lll}  were  compared  (Table  2). 


TABLE  2 

500  eV  ENSEMBLE  RESULTS 


Number  of  trajectories 

300 

150  (x's) 

150  (+'s 

#  R h  atoms  sputtered 

902 

457 

445 

AS  I 

3.  01 

3.  05 

2.  97 

Average  energy  of 

sputtered  Rh  (eV) 

13. 272 

13. 192 

13. 355 

1st  layer  yield  ratio 

0.  970 

0.  969 

0.  971 

Note:  the  x's  and  +'s  denote  the  separate  groups  of  impact 
pointe  delineated  in  Figure  5. 

In  addition,  the  angular  and  energy  distributions  of  the 
150  trajectory  simulations  were  identical  in  form,  and 
nearly  so  numerically.  Since  all  the  measurable  results 
were  very  similar,  it  was  deemed  prudent  to  conduct  all 
subsequent  simulations  using  only  150  trajectories  (those 
designated  by  an  x  in  Fig.  5). 

4.  Adatom  and  Vacancy  Placement 

Adatoms  are  placed  at  the  3-fold  bridge  position 
(Fig.  6)  of  a  clean  Rh(lll}  surface  by  designating,  within 
the  simulation,  their  relative  position  from  a  surface 
layer  atom,  as  measured  from  atom  center  to  atom  center. 
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This  corresponds  to  2. 187A  in  the  negative  y  direction  and 
1.56QA  in  the  positive  z  direction  from  the  surface  atom. 
For  simplicity,  an  adatom  will  be  designated  by  the  number 
of  the  surface  atom  with  which  it  is  associated  (e.g., 
adatom  108).  See  Figure  6  for  an  example  of  this 

designating  system,  and  Figures  3  and  4  for  the  atom 


numbering  system. 

Vacancies  ar 

e  treated  in  a 

m 

uch  simpler  fashion. 

At  those 

positions 

where  a  vacancy 

is 

desired, 

the  atom  is 

simply  removed  from 

the  lattice. 

The 

vacancy 

location  is 

designated 

by  the 

original  atom 

*  3 

number  ( e. 

g. .  vacancy 

52). 


A  perfect 
these  simulations, 
distortion  from  the 
lattice  due  to ’ the 
For  example,  the 
toward  a  vacancy. 


C.  POTENTIAL  FUNCTIONS 
1 .  Background 

At  present,  the 


potential  is 

relatively  good 

(  i.  e.  ,  in 

the 

region  of 

greater  ) ,  due 

•  to 

the  large 

available. 

Also, 

the  theory 

separations 

where 

nuclear 

always  assumed  in 
made  to  account  for  any 
position  of  atoms  in  the 
of  adatoms  or  vacancies, 
surrounding  atoms  inward 


knowledge  about  interatomic 
at  large  separation  distances 
equilibrium  separation  and 
amount  of  experimental  data 
is  satisfactory  at  very  small 
scattering  is  predominant. 


crystal  lattice  is 
Mo  attempt  is 
equilibrium 
introduction 
relaxing  of 
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a  serious  lack  of  information  in  the 


However,  there  is 
region  between  the  two  extremes,  which  is 

compensated  for  by  extrapolating  either  high  e 
equibrium  potentials  based  on  a  reasonable  physica 
The  potentials  used  in  sputtering  models  lie 
intermediate  region  and  are  therefore  very  approx im 

A  typical  two-body  interatomic  potential 
of  an  attractive  part  at  large  separation  appro 
minimum  in  the  region  of  equilibrium  separati 
becoming  repulsive  and  increasing  rapidly 
interatomic  separation  decreases  further  and  th 
shells  of  electrons  of  the  atoms  begin  to  overl 
exact  form  of  the  interaction  in  the  low  energy  reg 
differ  for  two  atoms  interacting  in  a  crystal 
diatomic  molecule  since  in  the  former  case  it  is  in 
by  atoms  in  its  neighborhood. 

Like  the  analytic  theories,  molecular 
simulations  treat  atom-atom  repulsive  interaction 
body  central  forces.  The  assumption  is  made 

forces  and  potential  energies  are  functions  only 
properties  of  the  atoms  or  ions  and  the  inte 

separation.  During  an  MI  timestep,  several  atoms 
exerting  forces  on  a  single  particle;  so  the  resu 
calculated  as  a  sum  of  the  individual  p 

interactions,  improving  the  calculation. 


normally 
nergy  or 
1  model. 

in  this 
ate . 

cons is t s 
aching  a 
on,  then 
as  the 
e  closed 
ap.  The 
ion  will 
and  in  a 
f 1 uenced 

dynamics 
as  two- 
that  the 
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rnuclear 
may  be 
ltant  is 
air- wise 
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Attractive  forces  are  included  in  this  model,  even 
though  they  have  little  effect  on  the  particle  dynamics. 
They  do,  however,  dominate  the  processes  by  which  atoms 
escape  the  surface;  so  their  magnitude  strongly  influence 
the  atom  yield  from  a  single  trajectory. 

The  cohesive  energy  and  atom  separation  of  the 
material  determine  a  well  depth  and  "width*,  establishing  a 
reasonable  approximation  to  a  pair -potential  well,  but  an 
attractive  force  range  must  be  assumed.  Prior  to  this 
study,  test  trajectories  were  run  using  near est -neighbor , 
NN ( 1 ) ,  interactions,  next  nearest -neighbor  inter act ions, 
NN(2),  etc.,  out  to  MN<4).  In  the  fee  lattice,  differences 
beyond  NN ( 2 )  are  very  minor,  so  normally  a  potential 
function  is  truncated  between  NN<2)  and  NN<3).  The  well- 
depth  of  the  attractive  potential  function  is  adjusted  so 
that  the  calculated  potential  energy  of  an  atom  in  the 
center  of  the  target  equals  the  heat  of  atomization  of  the 
target  material. 

2.  Functional  Form 

Since  the  early  days  of  the  statistical  theories  it 
has  been  predicted  that  the  sputtering  yield  would  be 
relatively  insensitive  to  the  form  of  a  particular 
interaction  potential  function  [Refs.  27,  56].  The 
critical  assumption  in  any  sputtering  calculation  continues 
to  be,  however,  the  choice  of  the  potential  functions. 
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The  potential  functions  used  in  this  simulation  are 


empirical  in  nature.  They  do,  however,  provide  a  more 
realistic  view  of  the  atomic  interactions  within  a  cascade 
than  potentials  derived  exclusively  from  theoretical 
considerations  which  themselves  contain  approximations. 
The  empirical  potentials  are  based  on  simple  analytic 
expressions  which  may  or  may  not  be  justifiable  from 
theory,  and  which  contain  one  or  more  parameters  adjusted 
to  the  model  or  deduced  from  crystal  equilibrium  data. 

In  general,  the  atom-atom  potentials  are  compound 
functions  consisting  of  a  repulsive  ”wall?t,  joined  smoothly 
to  an  attractive  section,  the  "well",  by  a  cubic  spline 
function.  Much  of  the  recent  work  has  been  done  with  the 
Born-Mayer  function  repulsive  potential  [Ref.  57], 

V (r )  =  A  exp  ( -Br ) 

The  parameters  can  be  thought  of  as  two  degrees  of  freedom 
in  the  function  corresponding  to  an  ion  size  and  hardness. 
Here  the  size  is  the  separation  at  1.  0  eV,  and  the  hardness 
is  the  pre-exponential  factor  A. 

Many  calculations  have  also  been  performed  with  a 
modified  Moliere  potential.  The  original  Moliere  form  is: 

V  ( r )  =  [  (Z1Z2e2/a)/ (r/a)  ]  j*(r/a) 

0(r/a)  =  [0.  35  exp(-0.  3  r/a)  +  0.  55  exp(-l.  2  r/a) 

+  0.  10  exp(-6.  0  r/a)  ] 

a  =  0.  8853a0/ [Z^ 1 /2  +  (Firsov  screening  length) 
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where  aQ  is  the  Bohr  radius,  and  and  are  the  atomic 

numbers  of  the  ion  and  atom  respectively.  Because  this 

function  is  intended  as  an  approximation  to  a  Thomas-Fermi 
potential  function,  Torrens  [Ref.  27]  has  shown 
justification  for  using  the  Firsov  screening  length,  a,  as 
an  adjustable  constant.  The  resulting  function  is  called  a 
modified  Moliere  function:  aJ  -  ka. 

The  attractive  portion  of  the  compound  potential 
function  has  the  Morse  potential  function  form: 

V(r)  =  De [exp<-2a(r  -  re)>  -  2exp<-  a(r  -  re)>) 

with  De  the  well  depth,  re  the  equilibrium  separation  of  an 
atom  pair  and  a  controlling  the  well  width. 

3.  Potential  Functions  Used  in  the  Model 

For  this  3tudy  the  atom-atom  potential  was  chosen 
a3  a  repulsive,  modified  Moliere  potential  joined  to  an 
attractive  Morse  potential.  The  composite  pair  potentials 
Vi -|  between  the  ith  and  j  th  atoms  separated  by  a  distance  R 
is  given  by: 

Vj^  j  =  A  exp  (-BR)  R  <  Ra 

vxj  r  c0  +  C i R  +  C2R2  +  C3R3  Ra  <  R  <  Rb 

vij  *  De [exp< -2a (R  -Re)>  -  2  exp<-a(R  -Re)>)  Rb  <  R  <  Rc 

vij  =  0  R  >  Rc 

The  functions  are  joined  smoothly  by  a  cubic  spline  at 
points  Ra  and  Rb.  The  modified  Moliere  function  is 
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truncated  at  Rc  so  only  NN ( 1 )  and  NN  ( 2 )  interactions  are 
inc  luded. 

A  Born-Mayer/Morse  composite  function  used  by 
Garrison  [Ref.  58]  is  provided  for  comparison  in  this 
study,  and  is  of  the  form: 

Vjj  =  A  exp  (-BR)  R  <  Ra 

Vij  =  C0  +  C4R  +  C2r2  +  C3R3  Ra  <  R  <  Rb 

Vij  =  De[exp<-2a(R  -Re)>  -  2  exp<-oc(R  -Re)>]  R^  <  R  <  Rc 

Vij  =  °  R  >  Rc 

The  interaction  of  Ar+  with  Rh  is  represented  m 
both  cases  by  a  purely  repulsive,  unmodified  Moliere 
function: 

vij  =  [ (Z1Z2e2/a)/(r/a) ]  *(r/a)  R  <  Ra 

vij  s  0  R  >  Ra 

The  parameters  that  define  the  potential  functions 
used  in  this  study  are  presented  in  Tables  3  and  4.  The 
modified  Moliere  potential  functions  for  both  RH74  and  RH7X 
are  created  by  multiplying  the  Firsov  length,  a,  by  0.740. 


41 


TABLE  3 

POTENTIAL  FUNCTION  PARAMETERS 


Rh-Rh 

Rh-Rh 

Rh-Rh 

Ar-Rh 

Designation 

RH74S 

RH7X7 

RHBMfl 

A  (  keV  ) 

0.  0 

0.  0 

54. 330 

0.  0 

B  (  A-  1  ) 

0.  0 

0.  0 

-5. 088 

0.  0 

D.  (eV) 

0. 8237 

0. 7595 

0. 8237 

0.  0 

R.  <  A) 

2.  750 

2.  750 

2.  750 

0.  0 

a  ( A-  1  ) 

1.  080 

1.  080 

1. 560 

0.  0 

R.  <  LU  ) 

0.  73 

0.  73 

0.  70 

2.  20 

Rb  ( LU ) 

0.  80 

0.  80 

0.  85 

2.  20 

Rc  (LU) 

2.  20 

2.  20 

2.  40 

2.  20 

k 

0.  740 

0.  740 

1.  0 

1.  0 

TABLE  4 

CUBIC 

SPLINE  PARAMETERS 

Rh-Rh 

Rh-Rh 

Rh-Rh 

RH7  4 

RH7X 

RHBM 

Co 

-4. 264x1 03 

-5.  613xl03 

4. 469xl03 

c, 

1 . 666x1 0* 

2.  198x10* 

-1. 548x10* 

c8 

-2. 159x10* 

-2. 859x10* 

1. 812x10* 

C3 

9. 299x 1 03 

1.  235x10* 

-7. 130xl03 

The  Rhodium-Rhodium  potential  functions  and  forces  of 
RH74  and  RH7X  are  compared  graphically  with  the  RHBM 
function  in  Figures  7-10.  The  force  between  rhodium  atoms 
is  determined  by  taking  the  derivative  of  the  analytic 
potential  functions.  The  two  force  functions  are  then 
joined  by  a  log-linear  spline,  which  causes  the  distinctive 
break  in  the  graphs  as  seen  in  Figures  8  and  10. 


Function  used  in  500  eV  Ar+  simulations. 

Function  used  in  the  2  keV  Ar*  simulations. 

Garrison's  Born -Mayer /Morse  potential  function 
[Ref.  581. 
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IV.  RESULTS 


A.  OVERVIEW 

The  results  presented  in  this  study  are  the  product  of 
a  series  of  QDYN  simulations  totalling  nearly  100  in 
number.  Two  basic  avenues  of  investigation  were  pursued  to 
determine  the  effects  of  small  scale  surface  defects  on 
sputtering.  They  include:  i)  the  effect  of  the  defect  size 
on  sputtering;  and  ii )  the  effect  of  the  distance  between  a 
single  surface  defect  and  the  representative  area.  The 
second  area  of  investigation  has  provided  a  measure  of  the 
surface  defect  density  required  to  induce  a  noticeable 
change  in  the  sputtering. 

The  variation  in  sputtering  as  a  function  of  defect 
size  was  examined  by  comparing  a  series  of  simulations  that 
included  a  clean  RhCill)  surface  and  surfaces  containing  1, 
2,  3,  and  7  vacancies  or  adatoms.  The  locations  of  the 
defects,  relative  to  the  impact  area,  are  depicted  in 
Figures  11  and  12.  Deviation  in  sputtering  was  similarly 
studied  by  comparing  the  results  from  a  clean  Rhflll) 
surface  with  simulations  that  included  a  single  adatom  or 
vacancy  at  various  distances  from  the  impact  area.  Both 
sequences  of  simulations  were  carried  out  at  bombarding  Ar+ 
energies  of  500  eV  and  2  keV.  The  investigation  of  several 
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side  issues  such  as  the 
and  containment  criteria 

Voluminous  amounts 
series  of  computer  sim 
portions  of  this  informa 
tractable  comparisons, 
vacancies  as  compared  to 
between  simulations  run  w 
keV.  The  comparisons 
experimentally  measurable 
angular  distribution 
discerning  trends  in  the 
discussed  to  account  for 
quantities. 

Table  5  provides  a  su 
used  in  the  majority 
conducted  for  this  study. 


comparison  of  potential  functions 
(see  Fig.  1)  were  also  conducted, 
of  data  were  collected  from  this 
ulations.  Presented  here  are 

tion,  assembled  into  a  series  of 
including:  i)  the  effect  of 

adatoms,  and  n)'  the  differences 
ith  ion  energies  of  500  eV  and  2 
are  made  amongst  various 

quantities  such  as  the  energy  and 
of  sputtered  atoms  to  aid  in 
data.  Various  mechanisms  are  also 
the  variations  observed  in  these 

mmary  of  the  operating  parameters 
of  the  Rh ( 1 1 1 } / Ar*  simulations 


TABLE  5 

QDYN  OPERATING  PARAMETERS 


Ion  energy  ( eV ) 

Ion  angle  of  incidence 
Lattice  dimensions 
Rh-Rh  potential  function9 
Ar-Rh  potential  function1 0 
Number  of  trajectories 
Representative  area 
pin  point  (atom  #) 


500 

90° 

18x5x12 

RH74 

Born-Mayer 

150 

51 


2000 

90° 

25x6x17 

RH7X 

Born-Mayer 

150 

108 


Modified  Mol i^re / Morse  composite  potential  function. 
See  Tables  3  and  4  for  function  parameters. 


1  0 


See  Table  3. 
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B.  SPUTTERING  YIELDS 


In  his  pursuit  of  containment  of  the  sputtering  process 
within  the  microcrystai  of  a  simulation,  Harrison  has  shown 
that:  "Only  the  first  few  layers  of  the  tar  got  are  needed 
t  o  com  pi  efei  y  os  pi  a  :i  n  flL  L  ox  per  i  men  tal  b  pu  1 1  er  i  n  q  res  u  Its , 
i  no  i  u  P  j  n .*  h  e  t  o  tal  yield" .  [Ref.  55]  Since  absolute 
containment  is  not  achieved  in  this  model  the  total  yield 
will  be  low  when  compared  to  experimental  results.  All 
other  computed  quantities,  however,  will  be  very  close  to 
their  Hinf inite  crystal”  values  for  the  lattice  sices  used 
in  these  simulations.  Relative  yields  are  therefore 
reported,  with  the  values  being  normalized  to  the  clean 
surface  results. 

Figures  13  -  16  detail  the  effect  of  a  single  defect  on 
the  normalized  yield  as  measured  by  the  defect's  distance 
from  the  center  of  the  representative  area.  The  center  of 
the  representative  area  was  chosen  as  the  reference  point 
rather  than  the  pin  point,  since  the  impact  points  are 
almost  uniformly  distributed  about  the  center. 

For  the  low  energy  case  a  significant  portion  of  the 
yield  is  initiated  by  the  recoiling  primary  knock-on  atom 
(PKA)  or  "target  atom”.  The  comparison  at  Ar*  energy  of 
500  eV  (Fig.  13)  indicates  that  a  vacancy  has  a  greater 
effect  on  the  yield  than  an  adatom  does.  This  is  to  be 
expected,  since  the  incoming  ion  has  the  opportunity  to 
transfer  its  energy  to  a  second  layer  atom  through  direct 
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collision.,  thereby  involving  atoms  much  deeper  in  the 
lattice.  Subsequent  collision  sequences  would  distribute 
the  energy  through  layers  deeper  in  the  microcrystal  than 
would  normally  be  expected .  The  likelihood  of  a  surface 
atom  having  enough  energy  imparted  to  it  to  eject  it  from 
the  surface  (about  25  eV )  is  thereby  decreased. 

The  presence  of  a  peak  at  1.5  LU  for  the  adatom  curve 
at  Ar+  energy  of  500  eV  (Figs.  13  and  lb),  corr esponding  to 
an  adatom  at  positions  42  and  6111,  cannot  be  reconciled  by 
a  similar  mechanism  and  further  study  of  the  detailed 
collision  sequence  is  required. 

At  Ar+  energy  of  2  keV  (Fig.  14),  significant  yield  is 
attained  from  both  the  recoiling  target  atom  and  the 
scattered  ion,  where  the  "ion"  yield  is  increased  by 
contribut ions  from  second  and  third  PKA's.  Again,  the 
vacancy  has  a  greater  effect  on  yield.  This  effect  is  most 
probably  caused  by  the  same  mechanism  described  above,  but 
in  this  case  the  energy  is  lost  from  the  microcrystal  by 
having  more  atoms  ejected  from  the  sides  and  bottom  of  the 
microcr yst al ,  as  compared  to  a  clean  surface. 
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tend  to  be  closer  to  the  surface 


than  an  ion  impacting  a 


clean  surface.  The  recoiling  atoms  would  then  have  an 
opportunity  to  transfer  a  greater  portion  of  their  momenta 
to  the  surface  atoms  directly,  causing  more  atoms  to  be 
ejected.  This  effect  is  not  present  at  500  eV  since  the 
Ar*  is  more  likely  to  be  reflected  off  the  adatom.  A  brief 
comparison  of  reflected  Ar+  is  presented  m  Table  6  to 
support  this  hypothesis. 


TABLE  6 

NUMBER  OF  REFLECTED  Ar  IONS 


Ion  energy  ( e V ) 

500 

2000 

Clean  surface 

58 

46 

1  adatom* 

96 

56 

adatoms  are  at 

the  same 

positions  relative  to 

iresentative  area. 

The  scattering  of  data  points  for  the  2  keV  adatom  case 
(Figs.  14  and  16)  was  investigated  in  some  detail.  The 
data  points  at  approximately  1.8  LU  correspond  to  adatoms 
at  positions  95  (normalized  yield  of  1.028)  and  120.  While 
these  adatoms  are  the  same  distance  from  the  center  of  the 
representative  area,  they  are  not  the  same  distance  from 
each  of  the  individual  impact  points.  Where  there  existed 
a  significant  difference  in  yield  per:  single  ion,  the 
sequence  of  collision  was  traced.  The  most  significant 
difference  between  the  two  trajectories  was  the  position  on 
the  surface  where  the  ion  actually  made  contact.  This  was 
caused  by  a  slight  difference  in  distance  between  the  ion's 
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starting  position  and  the  adatom.  The  long  range  of  the 
Ar-Rh  potential  function  (2.2  LU  for  cutoff)  caused  the 
incoming  ions  to  follow  different  paths  to  impact  which 
yielded  two  different  collision  sequences.  This  was 

confirmed  by  reducing  the  Ar-Rh  potential  function  cutoff 
to  1.6  LU,  at  which  time  the  trajectory  results  coincided 
exactly.  ‘Mo  trend  can  be  observed  for  the  dispersed  data 
points  as  a  whole  concerning  their  relative  positions  to 
one  another. 

In  general  the  graphs  asymptotically  approach  a 
normalized  yield  of  1.0.  This  corresponds  to  the  defect 
being  moved  farther  from  the  area  of  collision  events  until 
it  has  no  effect  on  the  trajectories  at  all.  The  distance 
at  which  the  curve  reaches  ±5  percent  of  the  clean  surface 
results  is  an  indication  of  the  density  of  defects  required 
to  visibly  affect  the  yield.  This  takes  into  account  an 
average  error  of  37.  as  measured  by  the  deviation  of  the 
curves  from  1.0,  and  uses  only  the  upper  data  points  of  the 
2  keV  adatom  curve.  Table  7  presents  the  approximate 

distances  as  determined  from  the  straight  line 

extrapolation  of  Figures  13  and  14,  and  the  surface  density 
corresponding  to  a  circle  enclosing  a  defect  at  that 
posi t ion . 
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TABLE  7 

SURFACE  DENSITY  OF  DEFECTS 
CAUSING  ±5%  DEVIATION  IN  YIELD 

Distance  of  defect  Surface  Density 

from  center  of  (defects/A  2  ) 

representative  area  ( LU > 


500  eV, 

1  vacancy 

1. 69 

0.  0153 

500  ev, 

1  adatom 

0.  57 

0.  0539 

2  keV, 

1  vacancy 

1 . 22 

0. 0237 

2  keV, 

1  adatom 

0.  63 

0. 0492 

To  add  some  credence  to  the  predicted  defect  density 
results,  simulations  were  run  at  2  keV  that  evenly 
distributed  the  defects  over  the  surface.  Figures  17  and 
18  depict  the  positions  of  the  vacancies  and  adatoms 
respectively.  The  closest  defect  to  the  center  of  the 
representative  area  is  actually  outside  the  distance  listed 
in  Table  7  in  order  to  provide  a  check  on  the  outer  limit 
of  results.  For  the  vacancies  a  raw  yield  of  957  was 
obtained,  versus  1022  for  a  clean  surface;  a  6.47.  decrease 
in  the  yield.  The  adatoms  produced  a  raw  yield  of  1001, 
which  is  only  a  2.17.  decrease  in  yield;  not  quite  the  57. 
that  was  desired. 

The  difficulty  with  the  adatom  case  is  that  the  next 
closest  position,  nearer  the  representative  area,  to  place 
an  adatom  would  have  located  it  within  the  representative 
area.  This  was  felt  to  be  too  biased  a  test  to  confirm  the 
results,  since  it  would  have  essentially  duplicated  the 
original  simulation,  and  was  therefore  not  conducted. 
Despite  the  lack  of  confirmation  of  the  adatom  results,  I 
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believe  the  the  defect  densities  necessary  to  produce  a  57. 


change  in  the  yield,  as  shown  in  Table  7, 
estimate. 


to  be  an  accurate 


C.  ATOMS  EJECTED  PER  SINGLE  ION 

With  the  advent  of  computer  modelling  the  ability  to 
determine  a  distribution  of  atoms  ejected  per  single  ion 
(ASI)  became  available.  This  quantity  cannot  be  determined 
experimentally  but  it  does  present  some  insight  into  the 
mechanics  of  the  system. 

Figures  19  and  20  present  a  series  of  comparisons 
within  a  single  ion  energy  regime  for:  a  clean  surface;  1 


and  3  vacancies;  a 

nd  1 

and 

3  adatoms. 

For 

the  500 

eV 

case. 

little  difference 

is 

seen  between 

the 

clean 

s 

ur f ace 

distribution  and 

both 

of 

the  adatom 

resu 

Its.  This 

is  as 

expected  from 

the  sputte 

ring 

yield  results. 

wher 

e  adatoms 

have  been  show 

n  to  have 

little 

impact  on  the 

total 

yield. 

Again,  as 

expected 

from 

the  sputtering 

yield 

results. 

vacancies  have  a  significant  effect  on  the  ASI 
distribution.  The  mass  of  the  distribution  tends  to  shift 
to  the  left  as  the  size  of  the  pit  formed  by  the  vacancies 
gets  larger,  with  the  low  yield  trajectories  becoming 
dominant . 


At  2  keV  much  the  same  behavior  can  be  seen  happening 
to  the  distributions,  with  the  same  effect  now  being 
observable  for  the  adatoms,  but  to  a  smaller  extent  when 
compared  with  the  vacancy  cases. 
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An  alternate  form  of  displaying  the  ASI  distribution 
which  accentuates  the  atom  ejection  process  is  shown  in 
Figure  21.  The  number  of  rings  in  each  circle  is 
proportional  to  the  ASI  value  for  an  impact  point  at  the 
center  of  the  circle.  The  resulting  pattern  indicates  the 
areas  of  the  representative  area  which  produce  high  yield, 
or  low,  by  the  relative  darkness  of  the  circle.  Since  the 
RHtlll)  surface  does  not  contain  any  true  channels,  there 
is  no  significant  difference  noted  in  any  atom  yield  per 
impact  point  diagrams  created  during  this  study,  Figure  21 
being  provided  as  an  example  only. 


D.  EJECTION  TINES 

Figures  22  and  23  provide  a  comparison  of  ejecti 
distr ibutions  for  Rh,  within  an  ion  energy  regime, 
are  affected  by  1  and  3  vacancies,  and  1  and  3  adato 
500  eV  there  is  a  slight  shift  to  lower  ejection  ti 
an  increasing  number  of  vacancies,  and  a  very  consi 
shift  for  the  adatom  cases.  The  shift  to  lower  e 
times  at  2  keV  is  barely  discernable  in  all  cases. 

The  most  significant  fact  that  can  be  seen  fro 
distributions  is  that  more  than  95%  of  the  atoms 
during  the  simulation  occur  prior  to  250  femptosecon 
<1  fs  =  10” 1 3  sec)  for  the  500  eV  cases  and  prior  to 


on  time 
as  they 
ms.  At 
mes  for 
der able 
jection 

m  these 
ejected 
ds 

300  f  s 


at  2  keV. 
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E.  LAYER  YIELD  RATIOS 

The  layer  yield  ratio  is  defined  as  Lrt  =  N* /N,  where  M* 
is  the  number  of  ejected  atoms  which  originally  resided  in 
the  nth  layer  of  the  target,  and  N  is  the  total  number  of 
atoms  ejected.  Lt  values  at  500  eV  and  2  keV  are  plotted 
in  Figures  24  and  25.  At  both  energies,  the  adatoms 
effectively  act  as  a  "new"  first  layer  causing  the  yield 
ratio  to  decrease  more  rapidly  than  the  vacancy  cases. 

For  all  cases,  except  at  low  energy  and  high  number  of 
defects,  the  ejected  atoms  come  primarily  from  the  first 
target  layer.  This  agrees  with  Harrisons  statement  that: 
"Under  all  conditions  85%  or  more  of  the  ejected  atoms  come 
from  the  surface  layer  of  the  target,  and  essentailly  all 
of  the  remainder  come  from  the  second  layer".  [Ref.  55] 
However,  as  the  number  of  defects  increase  in  the  low 
energy  case,  the  target  is  more  able  to  accomplish  a 
momenta  reversal  and  have  atoms  ejected  from  deeper  layers. 


F.  EJECTED  ATOM  ENERGY  DISTRIBUTIONS 

The  ejected  atom  energy  distributions  were  created  by 
tallying  the  energy  of  all  Rh  atoms  sputtered  from  the 
surface  of  the  target  (whether  the  Rh  was  single  atoms  or 
multimers)  into  energy  bins  0.5  eV  wide.  Figures  2S  -  28 
are  plots  of  the  individual  tallies.  An  attempt  to  plot  a 
smooth  curve  through  the  data  points  is  not  made,  since 
significant  noise  is  present  in  the  plotted  values. 
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Figures  26  and  27  present  information  obtained  from  2 
keV  simulations.  It  can  be  seen  that  there  is  little 
difference  between  the  energy  distr ibutions  as  the  number 
of  defects  increase.  These  distributions,  as  a  group,  can 
be  viewed  as  one  ensemble.  By  averaging  the  energy  values 
in  each  bin,  their  statistical  plausibility  could  be 
improved. 

Figure  28  was  compiled  from  simulations  run  at  500  eV 
with  300  trajectories  each,  since  those  completed  using  150 
trajectories  did  not  contain  enough  information  to  discern 
anything  beyond  the  noise. 

By  allowing  the  data  points  to  be  interpolated  by  eye, 
these  energy  distributions  can  be  seen  to  follow  a 


Maxwellian 

form, 

with 

a  maxima  at 

about  2. 

5  eV 

Additionally, 

there 

is  a  s 

light  secondary 

peak  at 

9.  3  e 

for  the 

2  keV 

cases 

and  one 

at  8. 0  eV  for 

the  500  eV 

graph 

At  this 

time 

no  explanation 

is  forwarded 

to  account  fo 

these  secondary  peaks. 

G.  EJECTED  ATOM  ANGULAR  DISTRIBUTIONS 

Since  a  major  thrust  of  Garrison's  work  [Ref.  58]  has 
been  to  explain  the  location  of  the  peak  of  the  polar  angle 
distributions,  as  compared  to  experimental  results  [Refs. 
48  -  50],  an  extensive  group  of  plots  (Figures  29  -  41)  has 
been  developed  to  compare  the  effect  of  ion  energy  and 
defect  size  on  the  distributions. 
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The  polar  angle  distributions  are  created  by  tallying 
only  those  ejected  Rh  atoms  with  energies  greater  than  5.0 
eV  and  summing  over  all  azimuthal  angles.  Discriminating 
against  the  very  low  energy  species  will  help  align  the 
results  more  closely  with  the  experimental  measurements. 

At  2  keV  the  shift  in  the  distribution  is  very  slight 
with  the  vacancies  forcing  the  distribution  to  smaller 
angles  and  the  adatoms  to  larger  ones.  This  effect  is  much 
more  pronounced  in  the  500  eV  results,  with  Figures  37,  40 
and  41  showing  the  shifts  most  clearly.  In  general,  as  the 
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H.  SPOT  PATTERNS 


An  alternate 

method  for  displaying 

angular 

diet  r 

ibut ion 

is  the  sputtered 

atom  spot  pattern. 

Unlike 

the 

angular 

distributions  di 

scussed  in  section 

G  above. 

this 

method 

reveals  both  the  azimuthal  and  polar  angle  dependence  of 
the  sputtered  atoms  simultaneously. 

Little  difference  was  noted  between  the  spot  patterns 
of  the  simulations  conducted  for  this  study,  therefore  only 
two  examples  are  presented  here  for  illustration  (Figure 


12  See  Figs.  11  and  12 
defects . 


for  relative  position  of 
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42).  Of  note  is  the  expected  six-fold  symmetry  of  the 

pattern  charact er istic  of  a  fcc(lll). 

These  graphs  were  produced  by  plotting  the  ejection 
angles  of  all  sputtered  atoms.  To  obtain  spot  patterns 
more  similar  to  those  experimentally  obtained,  such  as  SIMS 
data,  many  more  trajectories  w ould  be  required  (on  the 
order  of  1000).  Atoms  ejected  with  energy  less  than  5.0  eV 
would  then  be  eliminated  to  produce  a  sharper  image. 

I.  BORN-MAYER  POTENTIAL  FUNCTION 

The  majority  of  theoretical  work  on  the  R h  ’  1  i  1  }  'A r * 
system  has  been  conducted  at  Pennsylvania  State  University 
by  B.  J.  Garrison  [Ref.  58],  in  conjuntion  with  the 
experimental  studies  of  N.  Winograd,  et.  al.  [Refs.  48- 
50].  Due  to  the  close  collaboration  between  Professor 

Harrison  and  Professor  Garrison,  it  was  considered  de 
rigueur  to  compare  the  results  obtained  using  the  Born- 
Mayer  Rh-Rh  potential  function  (RHBM)  with  the  modified 
Molit-re  potential  function,  which  is  used  throughout  this 
study  (refer  to  Chapter  3,  section  D. 3  for  details  of  these 
functions ) . 

Simulations  were  run  on  clean  surfaces  only,  at  ion 
energies  of  500  eV  and  2  keV.  Significant  differences  in 
results  are  apparent.  Summaries  of  the  graphical 

comparisons  are  listed  below,  with  the  statements  being 
couched  as  to  the  effect  of  the  RHBM  function  when  compared 
to  the  Moli^re  function. 
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1)  Figure  43.  Atoms  sputtered  per  single  ion  (ASI): 
At  500  eV  a  noticeable  shift  to  lower  ASI  is 
observed,  with  only  a  minor  shift  appearing  at  2 
keV. 

2)  Figure  44.  Atom  ejection  time  distribution:  At 

both  energies  a  shift  to  lower  ejection  times  is 
evident,  with  the  times  required  for  95X  ejection 
being  reduced  by  60  is. 


3  ) 

Figure  45. 

Ejected 
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However, 

some  details  ca 

5) 


det  er  mi ned . 


There  is  a  slight  shift  to  lower 


energies 

for 

both 

500 

eV 

and  2  keV 

corr espo 

nding  redu 

ct  ion 

in 

the  maxima 

4  ) 

F igures 

46 

and 

47. 

Ejected 

distributions:  These  figures  present  the  most 
dramatic  differences  between  the  potential 
functions,  with  the  RHBM  results  being  very  sharply 
peaked  at  25°  for  both  energy  regimes,  compared  to 
broad  distributions  peaking  at  larger  angles  for 
the  Moliere  results. 

Figures  48  and  49.  Atom  yield  per  impact  point: 
These  figures  are  presented  for  general 
information.  Ho  significant  differences  were  noted. 
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Table  8  presents  the  raw  yields  for  each  of  the  four 
simulations  compared,  and  a  percent  change  in  the  yields, 
to  emphasize  the  magnitude  of  the  effect  the  choice  of 
potential  function  can  have. 


TABLE  8 


COMPARISON 

OF  RAW  YIELDS 

Ion  Energy  ( eV ) 

500 

2000 

Yield  ( RHBM ) 

376 

1032 

Yield  ( RH74 , RH7X ) 

457 

1022 

Percent  change 

-18 

*1 

57 


V 


CONCLUSIONS 


While  this  study  has  made  no  attempt  to  compare  the 
simulation  results  to  experimental  results  that  are 
currently  available,  it  has  demonstrated  that  the  presence 
of  defects  on  single  crystal  surfaces  can  produce 
measurable  differences  in  the  quantities  that  characterize 
sputtering.  By  investigating  yet  another  variable  in  the 
simulation  of  sputtering,  this  study  has  laid  to  rest  some 
of  the  criticisms  leveled  against  the  simulations  in 
general,  e. g. ,  a  simulation's  lack  of  accuracy  due  to  the 
incomplete  modelling  of  the  system.  In  the  future,  the 

placement  of  vacancies  and  adatoms  may  be  considered  when 
any  system  is  modelled  in  order  to  provide:  i)  a  more 

accurate  prediction  of  experimental  results,  ii)  possibly 
greater  correlation  with  the  existing  experimental  results, 
and  iii)  greater  insight  into  the  actual  sputtering 
mechanisms  as  simulations  more  closely  reflect  actual 
sputtering  events.  These  statements  are  predicated  upon 

the  assumption  that  future  experiments  will  be  conducted 
with  the  precision  evident  in  the  EARN  experiments 
conducted  by  N.  Winograd,  et.  al.  [Refs.  48  -  50]. 

The  correct  choice  of  a  modified  Holier e/Norse  potential 
over  a  Born -Mayer/Morse  potential  was  verified  when  the 
results  of  the  simulations  were  compared  to  preliminary 
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EARN  results  [Ref.  58  3.  Garrison's  simulations,  using  the 
RHBM  potential  function,  fail  to  correctly  predict  the 
peaks  occuring  in  the  angular  and  energy  distributions, 
while  the  simulation  results  of  this  study  are  a  very  good 
correlation  for  the  angular  distribution  only.  Correlation 
of  the  energy  distribution  is  left  for  later  studies. 

The  choice  of  potential  function  does,  however,  cloud 
the  issue  of  the  effect  of  the  surface  defects.  The  choice 
of  an  appropriate  potential  function  is  still  considered  to 
be  a  "black  art",  since  little  is  known  of  the  behavior  of 
potential  functions  in  general  at  the  low  energies  at  which 
the  sputtering  processes  operate.  The  choice  of  an 
appropriate  function  only  has  to  fit  a  few  parameters,  such 
as  the  surface  atomization  energy  of  the  crystal.  Beyond 
that,  the  person  creating  the  model's  environment  has  wide 
latitude. 


The  magnitude  of  the  eff 
choice  is  at  least  as  great 
reported  here,  so  it  can 
incorrect  potential  function 
do  not  believe  that  this 
presented  in  this  study, 
surfaces  used  in  experiments 
them  into  account. 

Future  studies  into  the 
simulations  should  increase  t 


ect  of  the  potential  function 
as  that  of  the  surface  defects 
be  argued  that  possibly  an 
was  chosen  fox'  this  model.  I 
is  born  out  by  the  results 
Surface  defects  exist  on  all 
and  simulations  should  take 


Rh { 111 } 

/ Ar+  system 

using  QDYN 

e  ion 

energy  to 

5  keV  to 
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currently  being  conducted  in  this 


match 

area. 

the 

obser 


the  experiments 
Also,  more  attention  should  be  paid  to  inves 
echamsms  that  may  be  responsible  for  the 
ed,  by  tracing  the  progress  of  the  collision 


t igat ing 
results 
events . 
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APPENDIX:  FIGURES 


This  appendix  includes  all  figures  referred  to 
body  of  this  report.  The  data  presented  in  these 
are  the  measurement  of  various  parameters  of  the 
rhodium  atoms.  The  energies  presented  on  each 
refer  to  the  impacting  argon  ion  energy  utilized 
simulation. 


in  the 
figures 
ejected 
graphic 
in  that 
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Fig.  1  Yield  versus  Lattice  Sice. 
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Fig.  3  18x5x12  Lattice  Numbering  System. 
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Fig.  4  21x6x17  Lattice  Numbering  System. 
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IMPACT  POINTS 
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Fig.  5  Impact  Point  Positions 
within  the  Representative  Area. 
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Fig.  &  Adatom  Position. 
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Fig.  7  Rhodium-Rhodium  Potential  Function. 
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Fig.  8  Rhodium-Rhodium  Forces. 
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Fig.  9  Rhodium-Rhodium  Potential  Function. 
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Fig.  10  Rhodium-Rhodium  Forces. 
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Vacancies 
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Fig.  11  Relative  Vacancy  Positions. 
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Fig.  12  Relative  Adatom  Positions 
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Fig.  13  Yield  versus  Distance  from 
Center  of  Representative  Area. 
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Fig.  14  Yield  versue  Dietance  from 
Center  of  Representative  Area. 
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Fig.  15  Yield  versus  Distance  from 
Center  of  Representative  Area. 
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Fig.  16  Yield  versus  Distance  from 
Center  of  Representative  Area. 
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Fig.  17  Distributed  Vacancy  Positions. 
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Fig.  IS  Distributed  Adatom  Positions. 
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Fig.  19  Atoms  per  Single  Ion  Distribution. 
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Fig.  20  Atoms  per  Single  Ion  Distribution. 
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Fig.  21  Atom  Yield  per  Impact  Point. 
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Fig.  22  Atom  Ejection  Time  Distribution. 
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Fig.  23  Atom  Ejection  Time  Distribution. 
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Fig.  24  Yield  Ratio  versus  Size  of  Defect. 
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Fig.  25  Yield  Ratio  versus  Size  of  Defect. 
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Fig.  26  Ejected  Atom  Energy  Distribution. 
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Fig.  27  Ejected  Atom  Energy  Distribution. 
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Fig.  28  Ejected  Atom  Energy  Distribution. 
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Fig.  29  Ejected  Atom  Angular  Distribution. 
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Fig.  30  Ejected  Atom  Angular  Distribution. 
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Fig.  31  Ejected  Atom  Angular  Distribution. 
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Fig.  32  Ejected  Atom  Angular  Distribution. 
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Fig.  33  Ejected  Atom  Angular  Distribution. 
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Fig.  34  Ejected  Atom  Angular  Distribution. 
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Fig.  35  Ejected  Atom  Angular  Distribution. 
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Fig.  36  Ejected  Atom  Angular  Distribution. 
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Fig.  37  Ejected  Atom  Angular  Distribution. 
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Fig.  38  Ejected  Atom  Angular  Distribution. 
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Fig.  39.  Ejected  Atom  Angular  Distribution. 
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Fig.  40  Ejected  Atom  Angular  Distribution. 
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Fig.  41  Ejected  Atom  Angular  Distribution. 
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Fig.  42  Spot  Patterns. 
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Fig.  43  Atoms  per  Single  Ion  Distribution. 
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Fig.  44  Atom  Ejection  Time  Distribution. 
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Fig.  46  Ejected  Atom  Angular  Distribution. 
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Fig.  47  Ejected  Atom  Angular  Distribution. 
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500EV  RII(111)/AR<111>  CLEAN  SURFACE  500EV  RII(U1)/AR<111>  CLEAN  SURFACE  (RIIBM) 

ATOM  YIELD  PER  IMPACT  POINT  ATOM  YIELD  PER  IMPACT  POINT 


Fig.  48  Atom  Yield  per  Impact  Point. 
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2KEV  RII(111)/AR<111>  CLEAN  SURFACE  2KEV  Rll(lil)/AR<111>  CLEAN  SURFACE  (RHBM) 

ATOM  YIELD  PER  IMPACT  POINT  ATOM  YIELD  PER  IMPACT  POINT 


Fig.  49  Atom  Yielc.  per  Impact  Point. 
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